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We extend the recently developed strong coupling, dimensionally reduced Polyakov-loop effec- 
tive theory from finite-temperature pure Yang-Mills to include heavy fermions and nonzero chem- 
ical potential by means of a hopping parameter expansion. Numerical simulation is employed to 
investigate the weakening of the deconfinement transition as a function of the quark mass. The 
tractability of the sign problem in this model is exploited to locate the critical surface in the 
(M/r, pL /r, T) space over the whole range of chemical potentials from zero up to infinity. 
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1. Introduction 

Solving QCD via lattice simulations has proven to be a formidable problem, even unsurmount- 
able, as of now, as soon as the quark chemical potential = /Ib/3 is switched on, i.e. if the finite- 
density part of the phase space is under study. Hence, various alternative approaches have been 
developed to gain knowledge; some are based on /i = and extrapolate to finite jj., while other 
rely on building effective models and approximate descriptions that capture, to a certain extent, the 
basic dynamics of the system under study. 

A dimensionally-reduced effective theory based on strong-coupling expansion was introduced 
for the pure gauge sector in [[T]] : it offers robust predictive power in locating the thermal transition, 
can be improved order by order in a systematic fashion and studied numerically with relatively 
small efforts. Here we report on the inclusion of heavy fermions in the theory, implemented through 
a hopping-parameter expansion, and on the introduction of a nonzero chemical potential, with a 
sign problem well under control even at large /x. 

This contribution offers a sketchy overview on the subject and focuses only on some of its 
features: for a more detailed discussion, we refer the interested reader to ||^. 

2. Effective theory 

The theory under study comes from applying simultaneously strong-coupling and hopping 
parameter expansions to the Wilson action (on a lattice with temporal extent aN-c = l/T, lattice 
spacing a, and gauge coupling j8); it is then suitable to investigate, with the advantages of a simpli- 
fied, dimensionally reduced model, the heavy quark region of the QCD phase space. 

It is possible to integrate out the spatial links by means of a strong coupling expansion, which 
results in an action whose terms are each given by an effective coupling (function of the original 
parameters jS , A^^, M ^i^d the hopping parameter k) and consist of Polyakov loops L, = TrW, = 
Trn^=i Uo{i,T). The partition function thus correctly reproduces the Z3 centre symmetry of the 
gauge sector as well as its breaking by the introduction of a finite quark mass M. In practical 
applications, we restrict ourselves to just a few terms in such an effective action. 

A remarkable aspect of this theory is that the definition of its partition function is not expressed 
with an action linear in the couplings; this is due to the possibility of performing a partial resum- 
mation among certain classes of graphs, which appears to improve convergence. Also, this calls 
for a careful definition of the suitable observables to characterise the phase structure. 

Moreover, in the pure gauge case it has been already observed that the resulting phase transi- 
tion, albeit remaining first-order, is much weaker than in the linear (i.e. un-resummed) case, thus 
resembling QCD more closely. The critical effective coupling Ao for the quarkless theory can be 
translated into a table of Pd^t) by means of strong-coupling mappings, obtaining results close to 
those of full 4D simulations, allowing for a continuum extrapolation, which yields Tc = 250(14) 
MeV (Fig. |T]). 

Heavy quarks enter the model through an expansion in the hopping parameter k; this results 
in a sum over closed loops which translates to a series expansion in mixed powers of K and u 
(the latter being the first non-trivial coefficient in the character expansion of the Yang-Mills action, 
M(j8) = j8/18 + . . .). Partial resummations within classes of similar terms lead to writing the quark 
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Figure 1: Continuum limit of the pure gauge transition temperature Tc , from the effective theory. 



contribution in the form of a determinant. If a quark chemical potential /i is turned on, each loop 
will pick up an additional factor e"^'^^ = e^l^ raised to the power of its winding number. 

When rewriting the model in terms of L,, a "potential" term appears, encoding the reduced 
Haar measure on the group and the Jacobian from expressing each L, as: 

ALe'^ =ddd^e^'' , L(0, 0) = e'^ + ^^'^ + , ^.2K(l) ^ 27 - 18|Lp + SReL^ - \Lf ; 

(2.1) 

the partition function studied has then the form: 



Z{X,h,h)= fUdL^e^^l n [l+2AReLiL*] ) ( ndet[(l+/iW,)2^'(l+/iW;)2^' 

-fx \ <ij> J \ ^ 



, (2.2) 



with effective couplings X{p,Nt,K) as given in [[I], §], and, to leading order, h = (Ixe"^)^^ and 
h = {iKe^"^)^^ (expressed to higher orders in [§]). The number of flavours Nf is from now on set 
to one, although thanks to the small values of h,h involved a linear approximation can be safely 
used to restore h — > Nfh. 

One can express the fermion part entirely as a function of L, L* : 

Q, = det[{\+hW,){\+hW])f = [i\+hL, + h^L; + h^){l+hL;+tL, + t)]^ . (2.3) 

Nonlinearities aside, the above partition function can be compared to a three-state Potts model 
with a spin-spin interaction (-^ A) and an external magnetic field (~ h,h) acting on each spin: from 
knowledge of the Potts case [Q, ||], that has the same symmetry pattern, we expect a phase structure 
in {h,X) at zero chemical potential (meaning h = h) as depicted in Fig. ^ (left). If the chemical 
potential is switched on, we have h^h, but the qualitative shape of the phase structure should not 
change. In the latter case, however, we use for convenience the "reduced" he-^l'^ = h. 

3. Numerical results 

The numerical investigation was performed with a Metropolis algorithm directly implement- 



ing the partition function Eq. 2.2: the sign problem appearing as soon as 7^ is treated in the 
usual way by folding the phase into the observable and updating according to the norm of the 
configuration weight. It turns out that on systems as large as = 24^ the average sign is well 
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Figure 2: Left: Expected phase space for zero chemical potential. Right: Pseudocritical line for ji =0. The 



line is a fit to Eq. 3.2 Also shown is the critical point. 



larger than zero up to values of /x/r of three or more, thus posing no big trouble.^ Each datapoint 
produced represents a statistics of about 10^ configurations, analysed with the binning technique in 
order to estimate uncertainties meaningfully. 

Besides the usual observables, suitable for an action linear in its couplings, we also use their 
"nonlinear" counterparts (apart from trivial factors, they reduce to the former for — > 0): 



(3.1) 

From these observables, the susceptibility and the Binder fourth cumulant have been built as 
Xo =^'^{{0^) ~ (O)'^) and B4 = -^^^q^^p-- The main goal of this work is to map the phase 
structure in the {h,X, j) space: first, the case of zero chemical potential is studied, then we intro- 
duce a real /x. 

3.1 Zero chemical potential 

The investigation proceeds in two steps: first, the pseudo-critical line Xpc{h) is mapped, sub- 
sequently its critical point {Xc,hc) is located. The pseudo-critical line is found by fixing six val- 
ues of 0.0002 < h < 0.0012, and for each value by performing a A-scan at various system vol- 
umes, identifying four volume-dependent pseudocriticality estimators (extrema of susceptibility 
and Binder cumulant of En^jQim)- Then, for each of those estimators, an infinite- volume extrap- 
olation Xpc{h,Ns) = Xpc{h) + ci{h)N^" gives a thermodynamic limit which we find to mutually 
agree. The whole pseudocritical line is parametrised as 

^pcih) = Xo-aih . (3.2) 

A fit to the six points works well (Fig. right) and gives ai = 1.797(18) and Ao = 0.18805(1), 
the latter roughly reproducing the pure-gauge critical point found in [|T]]. ^ The linearity can be 
explained by a first-order expansion, in the small couplings involved, of the free energy of the 
system [^]. 

'Also, the configurations yielding a minus sign from the gauge part are extremely rare and can be in fact ignored at 
these system volumes and in the region of parameter space of interest 

^The slight deviation between the two results is due to the small volumes used in M\ for the determination. 
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Figure 3: Effect of the rotation {E,Q) [E' , Q') in locating the jj. =0 critical line. Top panel: original 
observables. Bottom panel: rotated observables. On the left the (normalised, centred) histogram is shown 
for the sample case Ns = 20, h — 0.000742, on the right the behaviour of the third moment of Q (Q') is 
plotted for various system volumes. Vertical lines mark the critical he- Note that the largest volumes show, 
for the rotated Q', a third moment essentially zero around the critical point. 



In order to locate the critical point, we switch to the nonlinear observables, Eq. |3.1| , and mea- 
sure them along Xpc{h) at various system sizes up to Ns = 24. Close to the critical point, and 
focusing on Q, we expect the following scaling laws for the susceptibility and Binder cumulant: 

Xq = Nl'^fjc^ (x) , B4,e = fB,,Q {x) ■ x={h- h,)Ny^ , (3.3) 

with critical indices dictated by the three-dimensional Ising universality class, i.e. 7/v ~ 1.962, v ~ 
0.6302. Moreover, universality also implies fs^QiO) ~ 1.604. Writing f{x) as a series in x, the 
susceptibility and Binder cumulant data were fitted to the above expectation keeping A^^ > 20, 
with rather stable results against different truncations for f{x), different scaling windows, fixing or 
leaving free the critical indices, and we get the final values {Xc,hc) = ^0.1 8672(7), 0.00073 l(40)y 
Another, more rigorous method is employed in ^ to identify the critical point (see also |Q] 
for an application to a model similar to ours). The two-dimensional distribution of (£", Q) is subject 
to a rotation — )• {E',Q'), with zero covariance; then, the critical line is defined as the locus where, 
in the thermodynamic limit, the third moment of the centred marginal distribution of Q' vanishes, 
^ 3 = 0. We explicitly verified that, around the critical point and for large enough volumes, the 

rotated Q' essentially satisfies this requirement along the line identified as described above (Fig. |3|). 

3.2 Real chemical potential 

As already observed, with the sign problem well under control for our purposes, we basically 
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Figure 4: Left: slope ai of the pseudocritical line as a function of IJ./T, with its cosh(jU/r) description. 
Right: the very weak dependence of on the chemical potential. The point labelled "Ac(0), reweighting" 
was found with the same technique as for finite as a cross-check. The slight drift in Xc is explained by 
a finite-size correction whose amplitude decreases as fi grows: indeed, the more accurate determination 
("scaling") for zero chemical potential and the large-/i values agree very well. 



repeat the analysis performed for ;U = at several values of /i/r up to 3.0. The only difference is 
that, for each chemical potential, we generate data only at a single point and then reweight 

all results to a whole 2D grid of points (the reweighting factors are complicated by the nonlinear 
formulation Eq. [2.2[ , but if one knows the target couplings in advance the table of weights can be 
prepared as the configurations are explored by the Monte Carlo). 

With the same statistics as for n = 0, for each chemical potential the 2d grid of values 
64,2 (A, /j) was scanned for the line of local minima: the largest- volume result was taken as the 
pseudocritical line and fitted to Xpc{h;pL /T) = Xq{}Jl/T) — a\{pi /T)h; again, a linear relation was 



sufficient, and turned out to be a constant compatible with the one in Eq. g^. Furthermore, 
it can be argued that, neglecting higher-order corrections, the slope of the curve depends on jx/T 
as a\{pL/T) = Ccosh(/x/r), a behaviour that was confirmed numerically with C = 1.814(3), in 
agreement with the = slope (Fig. ^, left). 

As for the critical point determination, the pL =0 fits encouraged us to define hc{pi/T) as the 
value where B^ q = 1 .604, with an uncertainty estimated from the difference between this definition 
and the h at which the B^^ q for the volumes 22^ and 24^ cross each other. Remarkably, the critical 
Xcin/T) shows little or no dependence on the chemical potential (Fig. ^, right), which allows to 
rewrite the parametrisation of ai (/x/r) as 

cosh(/i/r) 

A fit of the measured points to the above curve works indeed well, giving D = 0.00075(1) in 
full agreement with the = result (Fig. ^, left). There is, however, a slight deviation from the 
above law, that we ascribe to higher terms of the expansion in (/i/r)^ which, to first order, led 
to parametrising ai{}x/T). Remarkably, a similar phenomenon occurs in the simpler case of an 
effective theory constructed with the three-state Potts model, as can be seen in Fig. ^, right. 

4. Conclusions and outlook 

With the knowledge of the curve h{iJ./T) we can use the heavy-quark approximate relation 
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Figure 5: Left: the critical curve hc{^/T) along with the best-fit to Eq. 3.4 for our effective model. Right: 
the equivalent figure in the case of the three-state Potts model 1^: the fit curve, restricted to jx/T < 0.5, 
shows that the same phenomenon of large-/x overestimation of the critical h occurs in both models as a 
(tiny) deviation from the assumption of constant Xc{pL/T). We plotted the data in terms of h for ease of 
comparison. 



h = Nfexp{—M/T) and obtain the critical surface in the upper-right corner of the Columbia plot; 
in particular, for the sake of comparison with existing literature [Q, we quote here the /i = 
values of Mc/T and K"c(A^t = 4) for A/y = 1,2,3 respectively: 

M 

^ = {7.22(5), 7.91(5), 8.32(5)} , k-,.(A^t = 4) = {0.0822(11), 0.0691(9), 0.0625(9)} . (4.1) 

By carrying on the expansions to higher orders, a more careful analysis of the feasibility of a con- 
tinuum limit can be performed; this program is indeed illustrated in [^, along with the study of the 
imaginary-/^ side of the phase space. Another interesting direction of this study is the investigation 
of the low-temperature, large density limit of the theory, besides, of course, the attempt to lower 
the fermion mass as much as possible, within the range of applicability of the hopping expansion. 
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